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Abstract 

We present the preliminary results of a search for leptoquark (LQ) pairs using 
(85.2 ± 3.7) pb _1 oi pp collider data collected by the D0 experiment at the 
Fermilab Tevatron from 1994-1996. We observe no evidence for leptoquark 
production and set a limit on a{pp —* LQLQ — » vv + jets) as a function of 
the mass of the leptoquark (Mlq), assuming BR(LQ — > z^g)=100%. At the 
95% confidence level, we exclude scalar leptoquarks for Mlq < 99 GeV/c 2 , 
and vector leptoquarks for Mlq < 178 GeV/c 2 . 
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The observed symmetry between the lepton (I) and quark (q) sectors suggests the ex- 
istence of a force connecting the two that is mediated by particles that couple directly to 
both leptons and quarks, and are therefore known as leptoquarks (LQ). Leptoquarks arise 
naturally as the vector bosons Jl| or Higgs particles || of a Grand Unified Theory [|IJ; as com- 
posite particles 0; as techniparticles ||; or as R-parity violating supersymmetric particles 

Leptoquarks carry both color and fractional electric charge. The Fermilab Tevatron can 
produce pairs of leptoquarks through the strong process pp — > g —>■ LQLQ + X with a 
production cross section that is independent of the coupling for scalar leptoquarks, but not 
for vector leptoquarks. In this study, we consider the specific cases of Yang-Mills coupling 
(YM), Minimal Coupling (MC), and the coupling resulting in the minimal cross section 

Decay between generations is theoretically possible; however, the limits from flavor- 
changing neutral currents imply that leptoquarks of low mass (O(Tev)) couple only within 
a single generation f7| . Decays of leptoquark pairs result in one of three possible final states: 
l^Pqq, ^uqq, and vvqq. This analysis || studies the vvqq final state, assuming BR(LQ — > 
uq)=100%. In a previous analysis of this final state 0, D0 set limits of M L q > 79 GeV/c 2 
for scalar leptoquarks, and M L q > 145 GeV/c 2 , 160 GeV/c 2 , and 205 GeV/c 2 , for vector 
leptoquarks with couplings that yield the minimum cross section ((J min ), Minimal Couplings 
(MC), and Yang-Mills (YM) couplings, respectively 0. The analysis presented here uses 10 
times more data than the previous analysis. The CDF collaboration has conducted a search 
for second and third generation leptoquarks with BK(LQ — > uq)=100% and set mass limits 
of 123 (148) GeV/c 2 for second (third) generation scalar leptoquarks and 171 (199) GeV/c 2 
and 222 (250) GeV/c 2 for second (third) generation vector leptoquarks with MC and YM 
couplings, respectively [|TTJ|] . 

The D0 detector ]TT|] consists of three major subsystems: An inner detector for tracking 
charged particles; a uranium-liquid argon calorimeter for measuring electromagnetic and 
hadronic showers, and a muon spectrometer. The jets measured with the calorimeter have 
an energy resolution of approximately a(E)=0.8\/E (E in GeV). We measure the missing 
transverse energy (]ftj>) by summing the calorimeter energy in the direction transverse to the 
beam. The measurement has a resolution of <r=1.08 GeV + 0.019(X|.&r|) (Et in GeV). 

We use an event sample defined by the selection criteria: 2 jets with Et > 50 GeV; $?> 
40 GeV; A0(jet, $ T ) > 30°; and AK(jet, jet) > 1.5, where ATI = ^(Ar/) 2 + (A0) 2 , r] is 
the jet pseudo-rapidity, and is the jet azimuthal angle. These criteria select events with 
high trigger efficiency Backgrounds arising from W or Z boson production are reduced by 
rejecting events with isolated muons or highly electromagnetic jets. We reduce cosmic ray 
backgrounds by rejecting events with jets containing little electromagnetic activity. The 
integrated luminosity after removing events corrupted by accelerator noise and detector 
malfunctions is 85.2 ± 3.7 pb _1 . 

The backgrounds in the sample consist of events with jets produced in association with 
a W or a Z boson, and events from top quark and multijet production. We use Monte Carlo 
generators to simulate the topologies of events with W or Z bosons or top quarks, and a 
GEANT simulation of the detector to predict the acceptance of these events. 

The W and Z backgrounds consist of processes containing only neutrinos and jets 
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(W — > ThV + jet, Z — > uu + 2 jets), processes with unobserved charged leptons 
(W — > l^v + 2 jets, Z — > jijj, + 2 jets, Z — > r h T[ + jet), and processes in which 
an electron is misidentified as a jet (W — > ev + jet, W — > r e u + jet). We use the 
PYTHIA Monte Carlo generator to predict the acceptances of the W/Z + jet processes, 



and the VECBOS Monte Carlo generator [|13|] to predict the acceptances of the W/Z + 2 
jets processes. We scale the generator cross sections to match the cross sections measured 
using the W and Z electronic decays. 

The top background consists of tt, tb, and tb production, where the top quark decays 
to an unobserved charged lepton, a neutrino, and a jet. We use the D0 measured cross 
section for tt production |TJj] and the calculated next-to-leading order cross section for the 



other processes fl5fl . We use the HERWIG generator |L6| to predict the acceptance of the tt 



process, and CompHEP |T7fl to predict the acceptances of the tb and tb processes. 



The multijet background arises primarily from 2 sources: Vertex mismeasurement and 
jet energy loss. To reduce the number of events with mismeasured vertices, we use the 
central drift chamber (CDC) to associate charged tracks with each central high p? jet. The 
tracks are used to determine the jet vertex position which is required to be no further than 
15 cm from the event vertex position. The latter is determined using all of the tracks in the 
event. We reduce the number of events with significant jet energy loss by requiring that the 
angle between the 1$ T and the jet with the second highest measured be greater than 60°. 

To predict the multijet background remaining in our sample, we use the sample of 
events whose jet vertex position deviates by 15 cm to 50 cm from the event vertex po- 
sition. We normalize this sample to the event sample using a multijet dominated sample 
(A0(jet 2, $ T ) < 60°). The upper bound of 50 cm provides the best agreement between 
the background prediction and the data for events with $ T between 30 GeV and 40 GeV, 
which is dominated by multijet events (Table |). Changing this value to 100 cm increases 
the multijet background prediction by 22% in this region, which we take as an estimate of 
the systematic error of the method. Table [TI| shows the total expected background and the 
observed number of events for the 2 jets + $ T data sample. 

To model the characteristics of leptoquark production, we use scalar leptoquark events 
generated with PYTHIA and vector leptoquark events generated with CompHEP. The cross 



sections for scalar leptoquark production have been calculated to next-to-leading order |L9 
while those for vector leptoquark production have been calculated to leading order pO 
The calculations use a renormalization and factorization scale of [i=Mlq, with theoretical 
uncertainties estimated by changing the scale to fi=M L g/2 and fi=2M L q. We use the lower 
cross section (/i=2M^q) in our optimization and in determining our mass limits. 

The analysis is optimized for the production of 100 GeV/c 2 scalar leptoquarks and 200 
GeV/c 2 vector leptoquarks, since leptoquarks with either of these masses would give a ~2cr 
excess if they exist. We use the JETNET neural network program, with the $ T and 
A</>(jet, jet) distributions as inputs for scalar leptoquarks, and the $ T and second jet pt 
distributions as inputs for vector leptoquarks. We show the neural network outputs and the 
chosen cuts for both of these masses in Fig. [I]. The cuts are chosen to maximize the inverse 
of the fractional error in the signal: 

Ni„ 

rirr ~ 



y/N lq +N background +AN? q +ANj* ai ' 



background 



6 



where Ni q and Nb ac k groun d are the number of signal and background events, respectively, and 
ANig and AN background are the associated uncertainties. We show the numbers of events 
after these cuts in Table |ITT[ 



Event Sample 


Number of Events 


Multijet 


162.8 ± 23.7 


W, Z, and top 


51.9 ± 7.0 


Total background 


214.7 ± 24.7 


Data 


224 


TABLE I. The expected and observed numbers of events 
$ T between 30 GeV and 40 GeV. 


in the multijet dominated sample of 


Background 


Number of Events 


Mnltiiet 


58 8 ± 14 1 ± 12 9 


(W — > ev) 4- iet 


51 9 ± 7 + so 7 


(W — ► tv) + jet 


46.3 ± 5.0 


(Z — > z^z/) + 2 jets 


36.1 ± 7.7 


(IT -> ^) + 2 jets 


18.7 ± 3.5 ±f;? 


ti^l^v + 4 jets 


10.6 ± 2.0 ± 2.3 


(IT -> ez/) + 2 jets 


8.3 ± 2.5 ±H 


(IT -> tv) + 2 jets 


5.6 ± 1.7 ±U 


t& -> / ± z^ + 2 jets 


2.0 ± 0.3 ± 0.2 


(Z — > tt) + jet 


2.0 ± 0.4 ±§J 


(Z — > ,u/x) + 2 jets 


1.7 ± 0.4 ±8i 


Total background 


242.0 ± 18.9 


Data 


231 



TABLE II. The expected and observed numbers of events in the 2 jets + $t sample. 



After applying the optimal cuts, we find that the observed number of events is consis- 
tent with the expected background, and that, consequently, we have found no evidence for 
leptoquark production. This null result yields the 95% confidence level cross section limit, 
as a function of leptoquark mass, shown in Fig. 0. We calculate the limit using a Bayesian 
method with a flat prior for the signal and Gaussian priors for background and acceptance 
uncertainties. The corresponding mass limits are 99 GeV/c 2 for scalar leptoquarks, and 178 
GeV/c 2 , 222 GeV/c 2 , and 282 GeV/c 2 for vector leptoquarks with couplings corresponding 
to the minimum cross section a m i n , Minimal Coupling, and Yang- Mills coupling, respec- 
tively. We summarize the various D0 mass limits as a function of BR(LQ — > V^q) for first 



7 



Z 60 
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PRELIMINARY 



Solid! Background 
Dashed: 100 GeV SLQ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Neural Network Output 



PRELIMINARY 



Solid! Background 
Dashed: 200 GeVVLQ 



1=L 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Neural Network Output 



FIG. 1. The neural network output for the data, for background (solid), and for leptoquarks 
(dashed). We show the optimization for 100 GeV/c 2 scalar leptoquarks (left) and 200 GeV/c 2 
vector leptoquarks with Minimal Coupling (right). We remove events to the left of the arrows. 



Leptoquark 


Ndata 


■^background 


n a 


a 95% (pb) 


100 GeV/c 2 Scalar 


58 


56.0 +& 1 


+2.1 


10.8 


200 GeV/c 2 Vector (MC) 


10 


13.3 t 2 2 -.% 


+2.6 


0.60 



TABLE III. The data, the expected background, the number of a excess that would be observed 
in the presence of a signal, and the 95% confidence level cross section limit. 
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generation scalar leptoquarks in Fig. [3] |J and for second generation MC and YM vector 
leptoquarks [^1| in Fig. |j. We note that the gap at low values of BK(LQ — > l^q) has been 
closed by this analysis. 




70 80 90 100 110 120 130 140 120 140 160 180 200 220 240 260 280 300 

M SUI (GeV/c 2 ) M M (GeV/c 2 ) 

FIG. 2. The 95% confidence level cross section limits as a function of leptoquark mass. We 
show the mass limits for scalar (left) and vector (right) leptoquarks. 
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50 100 150 200 250 
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FIG. 3. The D0 excluded mass vs. BR(LQ — ► eq) region for first generation scalar leptoquarks. 
The dark region is excluded by this analysis. 
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FIG. 4. The D0 excluded mass vs. BR(LQ — > fj,q) region for second generation MC (left) and 
YM (right) vector leptoquarks. The dark regions are excluded by this analysis. 
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